Interface enhancement of ferroelectricity in CaTiOa/BaTiOs superlattices 
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We carry out first-principles calculations for CaTiOs/BaTiOa superlattices with epitaxial strain 
corresponding to growth on a SrTiOa substrate, and consider octahedral rotations as well as ferro- 
electric distortions. The calculations are done as a function of electric displacement field, and both 
a macroscopic and a local electrostatic analysis are carried out. We find that strong octahedral 
rotations occur for TiOg octahedra sandwiched between CaO layers on both sides, but are strongly 
suppressed if either neighboring layer is a BaO layer. Due to the resulting enhancement of the ferro- 
electric instability in the BaO-neighboring octahedra, we find that overall the ferroelectric instability 
of the superlattice is enhanced by the interface. Thus, short-period superlattices in this system have 
a larger ferroelectric polarization than longer-period ones of the same average composition, contrary 
to the expected trend. 

PACS numbers: 77.22.-d, 77.22. Ej, 77.80.-e, 77.84.Lf 



Perovskite oxide superlattices are currently the subject 
of intense scientific interest Through variation of 

the choice of constituent materials and their order of de- 
position, as well as the epitaxial strain state via selection 
of the substrate, a wide variety of functional properties 
can be achieved. While the properties can be understood 
partly as the result of the behaviors of the constituent 
materials under the electrical and mechanical boundary 
conditions characteristic of the superlattice 0, in many 
cases the atomic and electronic features of the interfaces 
have been shown to play a critical role As a result, 
materials design via interface engineering is becoming an 
essential strategy in guiding experimental exploration of 
this enormous class of materials, and there is clearly a 
pressing need for improved atomic-scale understanding 
of the role played by the interfaces in determining the 
functional properties of perovskite superlattices. 

First-principles methods have been extremely useful 
in providing both qualitative and quantitative character- 
ization of interface effects in superlattices 0, @] . In most 
perovskite superlattices combining ferroelectric (FE) and 
paraelectric (PE) layers @, S i, m-Q , the FE instability 
and spontaneous polarization are found to be suppressed 
with increasing interface density for a given overall com- 
position. In an analysis based on layer polarizations de- 
fined in terms of localized Wannier orbitals [6| and using 
electric displacement D as the fundamental field variable, 
the properties of interfaces were found to be determined 
mainly by the adjacent atomic layers [5|,ll0|. In this local 
analysis, the suppression of ferroelectricity by interfaces 
was associated with the formation of an interface dipole 
pointing antiparallel to the local polarization. 

An additional feature of many perovskite oxides is a 
tendency for rotations and tilts to occur in the corner- 
sharing network of oxygen octahedra. For example, the 
crystal structure of bulk CaTi03 (CTO) is generated by 



large oxygen octahedral rotations; its nonpolar charac- 
ter is thus understood as the result of the competition 
between FE and antiferrodistortive (AFD) orders, which 



suppresses the FE instability in the AFD state [111 ]. In 
constrast, BaTi03 (BTO) is highly resistant to oxygen 
octahedral rotations and exhibits a robust FE state at 
room temperature [12]. Previous first-principles stud- 
ies of the CTO /BTO superlattice system have neglected 
the oxygen octahedral rotations 0, 0, EH, yielding an 
artificially strengthened ferroelectricity and an enhanced 
ground-state polarization 0, 0] which increases as the 
interface density decreases. It is quite surprising, then, 
that a slight increase in polarization from the 2:2 to the 
1:1 superlattice is observed in a recent experiment [lij ]. 

In this paper we show that an increased density of in- 
terfaces can in fact enhance the polarization of a superlat- 
tice, illuminating the experimental observation [13| |. The 
key to this unusual behavior is the role played by octahe- 
dral rotations, in addition to FE distortions, in determin- 
ing the energetics and polarization of the interfaces 
Working with superlattices comprised of CTO and BTO 
layers as our model system, we use first-principles cal- 
culations to demonstrate this novel effect and to clarify 
the reasons for its occurrence. We find that the octahe- 
dral rotations are large for TiO@ octahedra sandwiched 
between CaO layers, but are strongly suppressed if either 
neighboring layer is a BaO layer. The latter case results 
in a local enhancement of the FE instability and of the 
resulting polarization. As a result, short-period superlat- 
tices, in which the density of interfaces is higher, have a 
larger spontaneous polarization than longer-period ones 
of the same average composition. 

The details of the calculations are as follows. We carry 
out strained-bulk calculations on all distinct period-four- 
compatible superlattices built from CaTi03 (C) and 
BaTi0 3 (B) cell layers, namely BTO, 3B1C, 2B2C, 
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TABLE I: Computed tetragonal distortion (c/4ao) and spon- 
taneous polarization, in C/m 2 , of the CTO/BTO supercells, 
with (P,,) and without (P s ,o) octahedral rotations. 





BTO 


3B1C 


1B1C 


2B2C 


1B3C 


CTO 


Ps 


0.388 


0.333 


0.184 


0.161 


0.0 


0.0 


c/ia 


1.0595 


1.0359 


1.0114 


1.0110 


0.9939 


0.9798 


Ps,0 


0.388 


0.346 


0.266 


0.374 


0.447 


0.590 



1B1C, 1B3C, and CTO, stacked in the [001] direction. 
The in-plane lattice constant is fixed at ao=7.275 bohr, 
our computed equilibrium lattice constant for bulk cubic 
SrTiOa, corresponding to epitaxial growth on SrTiC>3. 
We use a 40-atom tetragonal supercell with lattice vec- 
tors of length \/2ao in the [110] and [110] directions 
and c « 4ao along the [001] direction. This allows us 
to consider both the FE distortion along [001] and an 
antiferrodistortive pattern of octahedral rotations about 
the [001] axis, resulting in P4bm space-group symme- 
try. We neglect the tilting of oxygen octahedra (rota- 
tions about an in-plane axis); we believe this is justi- 
fied because oxygen tilting requires a coherent pattern 
of tilts that would propagate into the BTO unit cells, 
where octahedral rotations are unfavorable. The struc- 
tural relaxations and electron minimizations were done 
at fixed electric displacement fields within the frame- 
work of density- functional theory in the local-density ap- 
proximation [15[ using the LAUTREC code package [lOj , 
which implements plane-wave calculations in the projec- 
tor augmented- wave framework fl6i | . We used a plane- 
wave cutoff energy of 80 Ry and a 4x4x1 Monkhorst- 
Pack k mesh. 

In Table I we present the spontaneous polarization P s 
and c/a ratio for each of the six superlattices consid- 
ered. As expected, P s increases with increasing BTO 
fraction. Comparing the two superlattices 2B2C and 
1B1C, which have the same overall composition, we 
find that the latter, with higher interface density, has 
a larger P s , contrary to the expectations based on FE- 
PE superlattices previously considered in the literature 
@, 0, 0] > but consistent with the experimental result re- 
ported in Ref. 14j. This enhancement in IB 1C, however, 
is not due to improper ferroelectricity as was the case for 
PbTi0 3 /SrTi0 3 described in Ref. @j. In both 1B1C and 
2B2C, the rotations suppress, rather than enhance, the 



polarization, as can be seen by comparing the polariza- 
tion in the presence of rotations with the polarization 
computed for rotations constrained to zero, given in the 
last row of Table I. Though the polarization in the ab- 
sence of rotations is much greater for 2B2C, reflecting 
the expected trend in which a higher interface density 
suppresses the polarization, the magnitude of the sup- 
pression by rotations is much greater for 2B2C than for 
1B1C, due to the large rotation angle in the CC layer. 
The net result is to reverse the trend, with the polariza- 
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FIG. 1: (Color online) Dependence of T1O2 layer polarizations 
(relative to the average of Ti02 planes in bulk BaTiOa and 
CaTiOs) on chemical environment in CTO/BTO supercells. 



tion in 1B1C with rotations being slightly larger than for 
2B2C. Application of a large electric field, as in Ref. [i~3j ]. 
is expected to favor states with higher polarization and 
smaller octahedral distortion, reflecting the competition 
between these two distortions; this could be investigated 
by application of the present approach but will not be 
further discussed here. 

To develop an atomic-scale understanding of the en- 
hancement of the polarization in 1B1C, we carry out a 
layer-by-layer analysis of structural and dielectric prop- 
erties of the superlattices @, @ . In previous work Q , we 
showed that the properties of each atomic layer depend 
only on the chemical identity of neighboring layers and 
on global parameters such as the in-plane lattice constant 
and the displacement field D (along z). We expect this 
"locality principle" to continue to apply when the AFD 
rotations are included. In Fig. [TJ we plot the layer po- 
larizations of Ti02 layers as a function of D for all the 
CTO/BTO supercells considered. For each D, we fully 
relax the structure, optimizing with respect to rotations 
of the octahedra in each layer. It can be seen that the 
layer polarizations are mainly determined by the chem- 
ical identity of the first-neighbor layers, just as was the 
case when octahedral rotations were constrained to zero 

We next consider the octahedral rotations and show 
that these, too, are functions of the local chemical en- 
vironment. In all of the CTO/BTO superlattices stud- 
ied here, we found that structures having rotations that 
change phase from layer to layer have lower energy than 
those with in-phase rotations. Henceforth we character- 
ize the structures by the absolute value 6j of the rotation 
angle for Ti02 layer j, with the understanding that the 
signs always alternate from layer to layer. 

The calculated dependence of the rotation angles on D 
is presented in Fig. [2] The results are grouped into three 
categories depending on whether the two first-neighbor 
AO layers are both CaO, one CaO and one BaO, or two 
BaO layers. It is clear that the octahedral rotations de- 
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FIG. 2: (Color online.) TiOe rotation angles as a function 
of D for TiOe octahedra sandwiched between (a) two CaO 
layers, (b) one BaO and one CaO layer, and (c) two BaO 
layers. 



pend strongly on the identity of the first-neighbor layers, 
while the influence of the second-neighbor layers is much 
weaker. A closer inspection shows that a large octahedral 
rotation of ^12-13° is observed whenever the TiC>6 layer 
is sandwiched by two CaO layers, while only a tiny rota- 
tion is found for the TiOg octahedra sandwiched between 
BaO layers. This is consistent with the properties of the 
parent materials, since rotations are strongly favored in 
bulk CTO while being suppressed in bulk BTO. We find 
that the D-field dependence of the octahedral rotations 
is generally weak in both cases. 

What is more intriguing, however, is the behavior of 
the interfacial layers of TiOg octahedra. These experi- 
ence a strongly broken mirror symmetry due to the pres- 
ence of BaO on one side and CaO on the other. The 
rotation angles are intermediate in size, ~4-6°, and are 
much more strongly dependent [l7j ] on D. (A similar 
effect due to second-neighbor asymmetry, barely visible 
for some of the curves in Fig. HJa) and (c), has a very 
much weaker amplitude.) It is this reduction in the an- 
gle that allows a larger polarization to develop in this 
layer, enhancing the overall polarization of the system. 

This behavior can also be analyzed by examining the 
tendency of a system towards ferroelectric behavior as 
measured by the inverse permittivity e _1 = dE/dD eval- 
uated at D—0 10], where E is the electric field. For a 
PE system e^ 1 (D = 0) is positive, and its approach to 
zero signals the proximity to a FE instability; for a FE 
system it gives a measure of the strength of the FE in- 



TABLE II: Reduced inverse dielectric permittivities e^T = 
eo/e computed at D=Q for the bulk and supercell structures 
studied. 





CaTiOa 


Rotations 


Rotations 


Superlatticc 


fraction 


allowed 


forbidden 


BTO 





-0.0191 


-0.0191 


3B1C 


0.25 


-0.0118 


-0.0108 


1B1C 


0.50 


-0.0030 


-0.0039 


2B2C 


0.50 


-0.0028 


-0.0074 


1B3C 


0.75 


0.0063 


-0.0070 


CTO 


1 


0.0157 


-0.0151 



stability that is independent of anharmonicity and thus 
complementary to the spontaneous polarization. 

In the third column of Table HI] we present the inverse 
permittivities e^ 1 {D = 0) for our bulk and superlatticc 
structures for the case that rotations are allowed. As ex- 
pected, the (strained) bulk BTO shows the most FE be- 
havior (strongly negative e _1 ), while the (strained) bulk 
CTO remains PE (e~ 1 > 0). The dielectric responses 
of the mixed superlattices are thus intermediate between 
those of bulk BTO or CTO, as expected. Consistent with 
the interface density dependence of the polarization, we 
find that the FE instability of the 1B1C superlattice is 
slightly greater than that of the 2B2C superlattice, and 
both of these are greater than the average of BTO and 
CTO, which corresponds to nBnC in the limit n — > oo 
and interface density equal to zero. 

The emergence of this interface enhancement effect is 
related to the coexistence of FE and AFD order in the 
system. To demonstrate this, the last column of Table 
fjll presents the inverse dielectric permittivities with oc- 
tahedral rotations constrained to zero. Here, we find 
that bulk CTO has nearly the same degree of FE in- 
stability as bulk BTO, while the presence of interfaces 
reduces the FE instability substantially for the 3B1C, 
2B2C, and 1B3C superlattices, and even more for the 
1B1C superlattice, with a two-fold increase in interface 
density. Comparing with the values when rotations are 
allowed, it is clear that the suppression of the octahe- 
dral rotation greatly increases the FE tendency in the 
supercells involving 50% or more CTO. 

From these comparisons, it follows that the competi- 
tion between AFD and FE orders is essential for the in- 
terface enhancement effect in this system. In fact, there 
is a local relationship between the octahedral rotations 
in a layer and their contribution to the FE instability. 
For the purposes of carrying out a local decomposition, 
we focus on the inverse of the capacitance per basal area, 
C" 1 = dV/dD, at D=0. Here V = LE is the poten- 
tial drop across the supercell, with L being the supercell 
height in the stacking direction and E the average elec- 
tric field. Apart from a small correction proportional to 
dL/dD, which we have verified has little effect on the 
analysis, this is equal to Le^ 1 , and thus can be used as 
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TABLE III: Inverse cell capacitance (Jm 2 /C 2 ) at D — for 
each TiCh-centered unit cell in the 2B2C supercell, for the 
case that rotations are allowed (P4bm symmetry) or forbidden 
(Pimm symmetry). 





r -\ 

W B(T)B 


°B(T)C 


C(T)C 


r -i 

°C(T)B 


°tot 


PAbm 


-0.416 


-0.361 


0.663 


-0.361 


-0.475 


P4mm 


-0.007 


-0.152 


-0.985 


-0.152 


-1.296 



a measure of ferroelectric instability similar to e _1 [Toj . 
We decompose C _1 locally as follows. First, for each AO 
or T1O2 layer in the superlattice, we define an inverse 
layer capacitance 



c- 1 - -fh +D^-^- 



(1) 



where hi and pj are the layer height |18| and layer polar- 
ization |6J] associated with each layer j. We then assign 
an inverse cell capacitance to the cell centered on T1O2 
layer i as 



-T1O2 



1 

2 C 



AO,i-l 



i c -l 

2 AO,i+f 



(2) 



Taking the 2B2C supercell as an example, we report 
the resulting inverse cell capacitances at D = in Ta- 
bic [Ell The last column of the table is C^l = Yli=i^i- 
For comparison, the results are given both for the case 
that rotations are allowed and when they are forbidden. 
As expected, the absence of the oxygen rotations results 
in a much stronger FE instability, as shown by the fact 
that C t ~J is much more negative than when rotations are 
allowed. The individual inverse cell capacitance show 
changes in Table I that are not always easy to interpret, 
because the inverse capacitance of a given TiC>2-centered 
unit cell does still depend fairly strongly on the identity 
of neighboring cells. For example, Cw T ^ B is not indepen- 
dent of environment, despite the fact that the rotations 
in that Ti02 layer are quite small. However, there is 
one huge change that stands out in Table I, namely, that 
Cc(T)C mverse capacitance of the cell containing the 
octahedra sandwiched between CaO layers, changes dras- 
tically when rotations are allowed. Here, the oxygen ro- 
tation angle is ~13°, similar to what it is in bulk CTO, 
and this induces an enormous suppression of the FE in- 
stability in that layer, switching the local tendency from 
a ferroelectric to a strongly PE one. 

These observations help explain the result presented 
earlier in the paper: that short-period superlattices can 
be more FE than long-period ones. We can expect that 
there are two competing influences: (i) the effect that 
typically occurs in FE-PE perovskite oxide superlattices 
in which the presence of an interface interrupts the FE 
order and reduces the FE tendency; and (ii) an effect 
that is driven by oxygen octahedral rotations and arises 



because of the severe suppression of ferroelectricity that 
occurs only for TiOe octahedra sandwiched between Ca 
layers on both sides. The density of such layers in an 
nBnC superlattice is (n — l)/2n, or l/2n less than in 
the average of bulk BTO and CTO. This l/2n reduction 
of a suppressing influence amounts to an enhancement of 
the FE tendency in short-period superlattices. Effects (i) 
and (ii) are of opposite sign but evidently have similar 
magnitudes, with (ii) dominating. This accounts for the 
increase in polarization for the increased interface density 
of the 1:1 superlattice relative to the 2:2 superlattice, as 
observed in a recent experimental investigation 14 1 . 

In conclusion, we have shown that oxygen octahedral 
rotations must be included for any realistic modeling of 
CTO/BTO superlattices. Using first- principles calcula- 
tions, we find that the octahedral rotations in a given 
layer depend strongly on the identity of the neighboring 
layers. In particular, the rotational angle in a Ti02 layer 
sandwiched between two CaO layers is found to be very 
large, while one between two BaO layers rotates hardly 
at all, similar to the tendencies of the corresponding bulk 
CTO and BTO parent materials. The competition be- 
tween rotations and ferroelectricity then plays out locally, 
and results in a novel behavior in which the ferroelectric- 
ity is found to increase systematically with decreasing 
n in nBnC superlattices, in contrast to what has been 
found when only the FE degrees of freedom are present. 
These results have critical implications for the design and 
application of short-period superlattice structures having 
enhanced or novel properties. 
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